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SUMMARY 

A  three  dimensional  finite  element  analysis  and  a  one  di¬ 
mensional  incompressible  flow  model  were  used  to  inves¬ 
tigate  the  interaction  between  the  internal  ballistics  and 
the  structural  grain  deformations  for  the  Sparrow  Mark 
52  rocket  motor.  The  results  confirm  earlier  investiga¬ 
tors’  conclusions  (based  on  two  dimensional  structural 
analysis)  that  a  mechanism  exists  for  pressure  spiking 
driven  by  choking  of  the  gas  fiow  within  the  bore  cavity 
due  to  grain  structural  deformations.  Significant  differ¬ 
ences  were  found  between  the  2D  and  3D  results,  how¬ 
ever;  the  3D  analysis  requires  a  higher  propellant  mod¬ 
ulus  to  prevent  choking  under  the  same  analysis  condi¬ 
tions. 

INTRODUCTION 

The  Sparrow  Mark  52  rocket  motor  had  a  history  of  oc¬ 
casional  “pressure  spikes”  during  the  ignition  transient 
in  test  firings  at  71  °C.  The  Mark  52  motor  (see  Figure 
1)  was  characterized  by  a  5-lobed  star  grain  configura¬ 
tion  and  a  relatively  constrained  flow  passage. 


Figure  1.  Motor  Geometry 


frequency  of  occurrence  of  pressure  spikes  in  Mark  52 
test  firings  was  increasing  with  motor  age,  so  we  were 
concerned  about  this  phenomenon  as  a  potential  age- 
out  mechanism  for  the  motor.  Since  accelerated  aging 
of  the  propellant  showed  a  tendency  to  soften  with  age, 
(see  Figure  2)  it  appeared  that  the  tendency  for  pressure 
spiking  should  be  age  dependent. 


Figure  2.  Propellant  Accelerated  Aging  Data 


The  purpose  of  the  work  reported  here  was  to  verify 
previously  reported  analyses  of  the  problem  with  an  im¬ 
portant  improvement:  while  the  earlier  analyses  used 
axisymmetric  (“2D”)  structural  models  to  represent  the 
propellant  grain,  we  used  the  TEXGAP-3D  program  to 
generate  a  detailed  three-dimensional  structural  model 
of  the  motor. 


Previous  analyses  of  the  Sparrow  Mark  52  motor  had 
indicated  that  the  anomalies  were  caused  by  an  inter¬ 
action  between  the  fiow  field  and  the  motor  structural 
response,  coupled  through  the  pressure/grain  deforma¬ 
tion  relationship.  The  nonlinear  interaction  between  the 
two  processes  can  become  unstable  in  these  motors  if 
the  propellant  modulus  is  low  enough,  causing  the  flow 
to  choke  near  the  end  of  the  grain  and  causing  the  pres¬ 
sure  spike.  At  the  time  this  analysis  was  conducted,  the 


ANALYSIS  PROCEDURE 

Figure  3  shows  the  overall  approach  that  we  used  to 
model  the  structural/ballistic  interaction  problem.  The 
key  parts  of  the  overall  model  are  the  geometry  model 
(composed  of  a  burnback  model  and  a  structural  defor¬ 
mation  model)  and  the  flow  model.  As  indicated  in  Fig¬ 
ure  3,  the  original  motor  geometry  is  altered  to  a  burn- 
back  configuration  by  the  burnback  model.  The  geom¬ 
etry  is  further  modified  by  the  structural  deformation 
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model  to  account  for  the  structural  deformations  (which 
are  defined  in  terms  of  a  deformation  amplitude  and  the 
motor  head  end  pressure).  The  final  modified  geometry 
and  head  end  pressure  are  then  fed  into  the  flow  model, 
which  computes  the  pressure  distribution.  The  struc¬ 
tural  deformation  model  within  the  geometry  model  is  a 
set  of  empirical  relationships  based  on  results  from  the 
TEXGAP-3D  analysis  and  other  structural  analyses. 


Figure  3.  Modeling  Approach 


Details  of  the  modeling  approach  are  discussed  later  in 

this  paper. 

ASSUMPTIONS 

The  major  assumptions  we  made  are  listed  separately 

below  for  the  two  main  parts  of  the  analysis:  the  struc¬ 
tural  analysis  and  the  flow  analysis. 

Structural  Analysis  Assumptions 

1.  All  materials  in  the  motor  behave  either  elastically  or 
in  a  quasi-elastic  manner,  characterized  by  effective 
viscoelastic  moduli  for  the  propellant  grain,  liner,  and 
insulation  and  by  constant  values  of  the  bulk  modulus 
of  compressibility  for  these  materials. 

2.  Linear  superposition  applies  for  all  structural  defor¬ 
mations;  furthermore  the  net  area  change  from  all 
contributions  (including  burning)  can  be  determined 
by  linear  superposition  of  the  components  of  area 
change. 

3.  The  structural  deformations  calculated  for  the  origi¬ 
nal  geometry  are  valid  for  the  burnback  geometry. 

Flow  Analysis  Assumptions 

1.  Friction  is  negligible. 

2.  The  gas  flow  is  steady  and  one-dimensional. 

3.  The  gas  can  be  treated  as  an  incompressible  fluid. 

4.  The  mass  generation  rate  per  unit  of  undeformed  sur¬ 
face  area  is  independent  of  pressure  and  location  with¬ 
in  the  propellant  grain. 


5.  Area  changes  due  to  structural  deformation  do  not 
change  the  mass  flow  rate  per  unit  of  undeformed  sur¬ 
face  area. 

STRUCTURAL  ANALYSIS 

We  used  a  number  of  different  structural  analysis  tools 
in  the  course  of  this  analysis.  The  structural  models 
used,  and  their  roles  in  the  overall  analysis,  are  discussed 
below. 


TEXGAP-3D  Model 

We  made  a  detailed  TEXGAP-3D  analysis  of  the  motor 
under  the  pressure  loading  shown  in  Figure  4. 


The  analysis  used  a  total  of  492  quadratic-displacement 
elements  to  model  a  basic-symmetry  section  (a  36-degree 
wedge)  to  capture  the  5-point  slotted  geometry  of  the 
motor.  The  elements  used  to  model  the  insulator,  liner, 
and  propellant  were  the  reformulated  elements  designed 
for  nearly  incompressible  materials.  In  addition  to  the 
internal  pressure  load,  the  boundary  conditions  included 
the  nozzle  ejection  load  and  the  igniter  ejection  load. 
The  deformed  geometry  of  the  aft  portion  (containing 
159  elements)  of  the  3D  model  is  shown  in  Figure  5. 


Figure  5.  Aft  End  of  TEXGAP-3D  Motor  Model 
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Figure  5  shows  a  primary  feature  of  the  grain  deflection 
of  the  Sparrow  Mark  52  motor  under  the  pressure  loading 
of  Figure  4;  the  propellant  grain  is  forced  against  the  aft 
dome,  causing  it  to  bulge  inward  and  constrict  the  cross 
sectional  area  available  for  gas  flow.  The  constriction  of 
the  flow  area  is  evident  in  Figure  6. 


Figure  6.  Deformation  of  Motor  Bore 


Reformulated  elements  were  used  for  the  propellant,  in¬ 
sulation,  and  liner.  Both  models  were  subjected  to  the 
same  boundary  conditions  as  the  3D  model. 

“QUICK-LOOK”  Code  Models 

The  Phillips  Laboratory  “QUICK-LOOK”  code  is  an 
upgraded  version  of  a  programmable  calculator  Lame- 
cylinder  structural  analysis  program  ^  which  was  ported 
to  the  Hewlett-Packard  HP-85  microcomputer.  The 
“QUICK-LOOK”  code  was  used  to  model  an  equivalent 
cylinder  for  the  slotted  grain  as  well  as  the  short  for¬ 
ward  circular  port  section  of  the  grain  which  surrounds 
the  igniter.  These  models  were  analyzed  under  general¬ 
ized  plane  strain  conditions  for  the  thermal  and  uniform 
pressure  loads. 

Material  Properties  Used 

The  material  properties  we  used  in  the  various  structural 
analyses  are  shown  in  Table  1.  The  properties  used  for 
the  nominal  analysis  (i.e.  propellant  modulus  of  1.55 
MPa)  are  essentially  the  same  as  those  used  by  Aerojet 
in  their  analysis  of  the  spiking  problem.  The  properties 
for  the  propellant,  insulation,  and  liner  were  considered 
to  represent  the  lower  limits  of  modulus  for  these  mate¬ 
rials.  The  lower  and  higher  propellant  modulus  values 
(0.689  MPa  and  3.45  MPa)  for  pressurization  were  in¬ 
cluded  to  allow  scaling  of  the  results  for  different  pro¬ 
pellant  modulus  values.  The  modulus  values  used  in 
the  “QUICK-LOOK”  analyses  were  based  on  effective 
modulus  calculations  for  a  specific  set  of  propellant  data 
discussed  later  in  this  paper. 


The  low  modulus  of  the  liner  material  (white  layer  in 
Figure  5),  permitting  a  large  shearing  deformation  of 
the  liner  layer,  is  a  key  factor  in  this  behavior. 


TEXGAP-2D/Approximate-3D  Models 
We  used  the  TEXGAP-2D/Approximate-3D  finite  el¬ 
ement  code  to  link  our  analysis  to  earlier  investiga¬ 
tions  which  used  two  dimensional  models.  Two  geom¬ 
etry  models  were  generated,  an  axisymmetric  equivalent 
cylinder  model  (basically  the  3D  geometry  with  all  ma¬ 
terial  between  the  slot  valleys  removed)  and  an  approxi¬ 
mate  3D  model  using  slot  elements  to  model  the  slotted 
region  of  the  grain.  Figure  7  shows  the  Approximate-3D 
finite  element  model. 


Complete  Grid  (Reduced  Size) 


PI  Approxlmate-3D  Elements 


Figure  7.  Approximate-3D  Finite  Element  Model 


These  models  used  153  and  237  elements  respectively, 
mostly  quadratic-displacement  isoparametric  elements. 


Table  1.  Mechanical  Properties  Used 
in  Structural  Analyses 


TEXGAP-3D 

TEXGAP-2D 

QUICK 

APPROX-3D 

LOOK 

PROPELLANT 

E,  MPa 

1.55 

0.689,1.55,3.45 

0.35, 1.55* 

/?*!* ,  MPa 

517 

517 

517 

LINER 

E,  MPa 

0.14 

0.14 

,  MPa 

6894 

6894 

INSULATION 

E,  MPa 

6.895 

6.895 

N/A 

,  MPa 

10300 

10300 

N/A 

CASE 

E,  MPa 

200000 

200000 

200000 

U 

0.32 

0.32 

0.32 

Used  for  thermal  deformation  only; 

...cd  a  =  3.6a;10“®cm/cm/°C.  AT  =  -6.6°C 

1  Poisson’s  ratio  (j/)  —  0.5  F/(6/3) 

where  P  is  the  Bulk  Modulus  of  Compressibility 
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Application  of  Structural  Analysis  Results 
The  primary  result  from  the  3D  model  was  the  radial  de¬ 
flection  at  the  bottom  of  the  slot  valley  under  the  pres¬ 
sure  load.  These  results  are  shown  in  Figure  8. 


Figure  8.  Radial  Deflection  at  Slot  Bottom 
from  TEXGAP-3D  Analysis 


In  addition,  the  deflections  at  nodal  points  on  the  bore 
surface  output  by  the  TEXGAP-3D  post  processor  were 
analyzed  using  a  specially  written  microcomputer  code 
to  determine  the  area  change  under  the  pressure  load. 

We  used  the  axisymmetric  equivalent  cyflnder  results 
only  for  comparison  with  TEXGAP-3D  and  the  Approx- 
imate-3D  model.  The  parameter  of  interest  is  the  radial 
deflection,  shown  in  Figure  9. 


Figure  9.  Comparison  of  2D  and  3D  Analysis  Results 


We  used  the  Approximate-3D  model  results  for  compari¬ 
son  with  TEXGAP-3D  and  also  (in  conjunction  with  the 
“QUICK-LOOK”  models  for  uniform  pressure)  to  de¬ 
velop  an  empirical  relationship  between  the  propellant 
modulus  and  the  radial  deflections.  As  shown  in  Fig¬ 
ure  9,  the  radial  deflection  determined  by  the  Approx- 
imate-3D  model  very  closely  matches  the  TEXGAP-3D 
results  in  the  aft  end  of  the  motor,  while  the  axisym¬ 
metric  model’s  radial  deflection  is  approximately  50% 


smaller  than  the  3D  model’s  deflection.  These  results 
emphasize  the  need  for  a  3D  grain  analysis. 

We  used  the  “QUICK-LOOK”  analysis  results  in  two 
ways:  (1)  to  define  the  grain  response  to  the  thermal 
load;  (2)  to  provide  the  uniform  pressure  component  of 
the  radial  deflection  at  the  slot  bottom  so  that  it  could 
be  subtracted  from  the  3D  deformations  under  the  to¬ 
tal  pressure  and  thus  determine  the  pressure  drop  com¬ 
ponent  of  the  this  deflection  (see  Figure  8).  The  same 
procedure  was  applied  to  the  peak  deflection  from  the 
Approximate-3D  model  to  determine  the  relationship  be¬ 
tween  the  propellant  effective  modulus  and  the  pressure 
drop  component  of  the  deflection.  The  resulting  empiri¬ 
cal  relationship  is 


(ArAp)max  =  [2.1234  -h  (3.684/E)](Ap/1.827)  (1) 

where  (ArAp)max  is  the  maximum  value  of  radial  de¬ 
formation,  in  mm,  at  the  slot  bottom  with  the  uniform 
pressure  contribution  subtracted  out,  E  is  the  effective 
elastic  modulus  of  the  propellant  in  MPa,  and  Ap  is  the 
pressure  drop  in  MPa  applied  to  the  grain.  This  empir¬ 
ical  equation  fits  the  Approximate-3D  results  (and  the 
single  3D  result)  with  a  maximum  error  of  1.2%. 

THE  GEOMETRY  MODEL 

The  geometry  model  modifies  the  motor  internal  geome¬ 
try  to  account  for  the  effects  of  propellant  burning,  struc¬ 
tural  deformation  under  thermal  load,  and  structural  de¬ 
formation  under  pressure  load.  The  cross  sectional  flow 
area  calculated  by  the  geometry  model  at  the  nominal 
analysis  conditions  {p^  =  12.07  MPa,  Ap  =  1.83  MPa, 
T  =  71  and  hum  distance  =  0.45  mm.)  is  shown 
in  Figure  10.  The  relative  magnitudes  of  the  various 
effects  are  evident  in  the  figure.  (The  alternative  calcu¬ 
lations  for  total  deformation- “3D”  and  “Rigid  Fin” -are 
explained  below. 


[Mstance  From  Head  End  of  Motor,  cm 

Figure  10.  Area  Variation  Along  Motor  Length 


Overall  Geometry  Calculations 

The  geometry  model  calculates  two  motor  geometry  pa¬ 
rameters  which  are  needed  by  the  flow  model.  These  pa¬ 
rameters  are  the  effective  radius  of  the  burning  perimeter 
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at  each  axial  station  and  the  cross-sectional  flow  area  at 
each  axial  station. 

The  burning  perimeter,  P,  is  represented  in  terms  of 
an  effective  radius,  a{z)  =  P/(27r).  The  effective  radius 
calculation  is  part  of  the  burnback  model.  The  cross- 
sectional  flow  area  is: 

A(z)  =  Aq(z)  AAb  +  AAp^  +  (2) 

where  A^^z)  is  the  original  flow  area,  AAb{z)  is  the 
change  in  flow  area  due  to  burnback,  and  the  remain¬ 
ing  three  terms  are  the  structural  deformation  compo¬ 
nents  due  to  thermal  load,  uniform  pressure,  and  pres¬ 
sure  drop,  respectively.  The  detailed  models  used  to 
generate  the  individual  geometry  change  components  are 
discussed  below. 

The  Burnback  Model 

The  burnback  model  consists  of  the  following  two  equa¬ 
tions: 

o(z)  =  00(2) +  1  (3) 

=  ’rj  [“(2)]  ^  -  [“0  (^)]  ^  }  (^) 

where  aQ^z)  is  the  original  burning  perimeter  and  x  is  the 
burn  distance,  which  is  specified  as  an  input  variable. 

The  Structural  Deformation  Model 

The  structural  deformation  model  computes  three  struc¬ 
tural  components  of  area  change:  the  thermal  load  com¬ 
ponent,  the  uniform  pressure  component,  and  the  pres¬ 
sure  drop  component. 

Thermal  Load  Component 

The  area  change  due  to  the  thermal  load  on  the  motor 
(temperature  soak)  is  estimated  from 

AAy  =  [(1  +  -  l]  -  2aATA^  (5) 

where  is  the  hoop  strain  in  an  equivalent-cylinder 
plane  strain  model  of  the  motor  cross  section,  is  the 
linear  coefficient  of  thermal  expansion  of  the  propellant 
and  AT  is  the  difference  between  the  soak  temperature 
and  the  motor’s  stress-free  temperature.  The  aAT  term 
arises  because  the  computed  strain  is  the  “stress-produ¬ 
cing  strain,”  i.e.,  ep  —  Ar/r^  —  a  AT.  The  strain  level 
was  found  to  be  nearly  independent  of  the  modulus  over 
the  range  considered. 

Uniform  Pressure  Component 

The  area  change  due  to  a  uniform  pressure  equal  to  the 
motor  head-end  pressure  is  calculated  from 

-A  (6) 

where  is  the  hoop  strain  in  an  equivalent-cylinder 
plane  strain  model  of  the  motor  cross  section  and  is 
defined  as  in  the  thermal  load  case  above.  Based  on  our 


plane  strain  analysis,  the  strain  level  is  nearly  indepen¬ 
dent  of  the  modulus  over  the  range  considered. 

Pressure  Drop  Component 

The  area  change  due  the  pressure  drop  through  the  mo¬ 
tor  bore  is  calculated  from  the  TEXGAP-3D  analysis 
results  and  the  uniform  pressure  plane  strain  analyses 
described  above,  the  specific  relationship  used  for  this 
calculation  is 

-  +  ^^Po  }  +  /(^)^’'Ap(^)] 

where  ArAP(^)  is  the  change  in  radial  displacement  at 
the  bottom  of  the  slot  valley  due  to  Ap{z)\  r^iz)  and 
Aq{z)  are  as  previously  defined;  and  AAp^  is  given  by 
Equation  (6).  The  empirical  function  f{z)  accounts  for 
the  effect  of  grain  fin  cross-section  deformations,  and  wa^ 
determined  from  the  deformed  and  undeformed  geome¬ 
tries  using  the  TEXGAP-3D  results  at  each  axial  station. 
As  used  in  the  structural  model,  ArAp(^)  was  calculated 
from 

where  Ar^^z)  is  the  TEXGAP-3D  radial  deformation 
with  the  uniform-pressure  component  Ar{z)\p=p^  sub¬ 
tracted  out,  and  Ar  is  the  input  strain  amplitude. 

Previous  analyses  of  the  Sparrow  Mark  52  pressure  spik¬ 
ing  problem  were  based  on  axisymmetric  (“2D”)  struc¬ 
tural  models  which  were  incapable  of  providing  the  three- 
dimensional  deformed  shape  of  the  propellant  grain.  As 
a  result  the  analysts  were  forced  to  rely  on  an  assumption 
termed  the  “rigid  fin”  assumption.  According  to  this 
assumption,  the  fin  cross-section  does  not  deform,  but 
simply  translates  toward  or  away  from  the  motor  center- 
line  to  conform  to  the  radial  deflection  of  the  equivalent- 
cyhnder  2D  model.  We  used  the  rigid  fin  assumption 
implicitly  in  our  thermal  load  and  uniform  pressure  mod¬ 
els,  where  we  expected  the  resulting  error  to  be  small  due 
to  the  essentially  plane  strain  boundary  conditions.  We 
considered  using  the  assumption  for  the  pressure-drop 
component  model  as  well,  since  it  would  have  simpli¬ 
fied  the  calculations.  To  determine  whether  a  rigid  fin 
assumption  was  acceptable,  a  “rigid  fin”  model  was  ob¬ 
tained  from  our  3D  model  by  setting  f(z)  =  0  in  Equa¬ 
tion  (10).  The  results  (see  Figure  10)  show  that  the 
“rigid  fin”  assumption  does  not  provide  an  accurate  area 
determination.  We  therefore  retained  the  empirical  f{z) 
representation  of  the  AA-to-Ar  relationship. 

THE  FLOW  MODEL 

The  theoretical  flow  model  we  used  was  originally  de¬ 
rived  by  Schapery  ^  and  was  applied  (in  a  somewhat 
simplified  version)  in  earlier  analyses  of  the  Sparrow 
Mark  52  pressure  spiking  problem. 

Based  on  the  assumptions  stated  earlier  along  with  the 
conservation  of  mass  and  momentum,  the  pressure  drop 
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equation  was  found  to  be 


Ap  =  Po  "  P  = 


(9) 


where 


C2=Ci 


(10) 


and 


As  before,  o  is  the  effective  radius  of  the  burning  perime¬ 
ter  and  Pq  is  the  motor  head  end  pressure.  The  other 
variables  are  L  =  motor  lengthy  At  =  nozzle  throat  area, 
p  =  gas  densitpy  ih  =  total  mass  flow  rate,  and  7  = 
ratio  of  specific  heats  for  the  gas,  the  gas  being  the  pro¬ 
pellant  combustion  products.  Note  that  if  we  assume  the 
variation  in  the  burning  perimeter  with  axial  distance  to 
be  negligible,  we  can  set  a  =  a  =  a^,  and  Equation  (9) 
simplifies  to 


Ap  = 


zdz 


(13) 


Equation  (13)  is  the  simplified  version  used  by  Schapery 
and  other  previous  analysts.  We  compared  the  results 
using  Equation  13  with  those  using  Equation  (9)  and 
found  them  to  be  in  close  agreement  except  that  the 
simplified  model  produces  an  abrupt  “dip”  in  the  pres¬ 
sure  at  the  sudden  expansion  in  the  forward  end  of  the 
grain,  while  the  pressure  calculated  by  the  more  complex 
model  varies  smoothly. 

We  implemented  Schapery’s  model  (Equations  (9)  to 
(12))  in  our  flow  model  by  assuming  a  piecewise-linear 
variation  with  z  for  both  a  and  A.  At  zero  burnback  dis¬ 
tance,  the  flow  analysis  results  show  a  reasonable  agree¬ 
ment  with  the  original  Aerojet  flow  analysis  results  (la¬ 
beled  “Assumed  Pressure”  in  Figure  11). 

By  repeating  the  flow  analysis  with  different  geome¬ 
tries,  we  were  able  to  show  that  the  difference  between 
the  Aerojet  results  and  our  zero-burnback  results  were 
caused  by  Aerojet’s  approximation  of  the  motor  port  ge¬ 
ometry  as  a  constant-area  duct. 

Figure  11  shows  that  the  expected  burnback  of  0.45  mm 
at  a  burn  time  of  0.36  ms  produces  a  significant  change  in 
the  pressure  distribution,  indicating  that  burnback  must 
be  included  in  a  realistic  analysis  of  the  pressure  spiking. 


Figure  11.  Initial  Flow  Analysis  Results 


Figure  12  shows  the  effects  of  area  changes  due  to  ther¬ 
mal  load,  uniform  pressure,  and  pressure  drop  on  the 
calculated  pressure  distribution  for  the  nominal  analysis 
conditions.  From  Figure  12,  it  is  evident  that  the  fi¬ 
nal  calculated  pressure  distribution  has  a  different  shape 
from  the  assumed  distribution,  which  violates  one  of  our 
original  assumptions.  While  we  did  not  directly  assess 
the  impact  of  this  lack  of  closure,  we  think  it  is  minor. 


Figure  12.  Flow  Analysis  Results 
at  Nominal  Analysis  Conditions 


RESULTS 

Taking  the  maximum  pressure  drop  from  Figure  12  as 
the  relevant  parameter  to  compare  to  the  nominal  pres¬ 
sure  drop  in  the  original  3D  structural  model,  we  were 
able  to  attempt  a  solution  of  the  nonlinear  problem  that 
defines  the  equilibrium  pressure  drop  for  the  motor.  The 
method  used  here  was  based  on  that  used  by  Click, 
Caveny,  and  Thurman  ^  to  evaluate  the  stability  of  a 
slotted-tube  propellant  grain  using  results  from  a  water- 
table  flow  simulation  and  a  propellant  grain  structural 
analysis.  The  method  used  was  to  first  plot  pressure 
drop  as  a  function  of  the  maximum  value  of  the  Ap  com¬ 
ponent  of  radial  displacement,  as  defined  by  the  flow 
analysis.  We  then  plotted  the  maximum  Ap  component 
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of  radial  displacement  versus  pressure  drop,  as  deter¬ 
mined  by  the  structural  model  (i.e.,  Equation  (1)),  and 
determined  where  the  curves  intersect.  Failure  of  the 
curves  to  intersect  corresponds  to  an  unstable  condition 
that  would  lead  to  very  large  deformations  and  choked 
flow  at  the  minimum- area  point.  This  condition  would 
produce  a  pressure  spike. 

Incorporation  of  Actual  Mechanical  Properties 
The  process  described  above  is  carried  out  in  Figure  13 
for  a  specific  set  of  propellant  modulus  values  based  on 
data  from  a  dissected  Mark  52  motor. 


Maximum  Inward  Radial  Deflection  (Pressure  Drop  Component),  mm 

Figure  13.  Analysis  Results  for  Ignition  at  71  ®  C 
(Head-End  Pressure  =  12.07  MPa) 


The  propellant  from  this  dissected  motor  was  very  soft. 
The  modulus  values  were  less  than  half  the  normally  ex¬ 
pected  values  over  the  entire  relaxation  spectrum.  How¬ 
ever,  since  this  propellant  (ANB-3109-2)  softens  with 
age,  the  modulus  data  from  this  motor  (SN  3600646)  was 
selected  to  obtain  a  conservative  evaluation  of  the  pres¬ 
sure  spiking  problem.  A  modified  power  law  relaxation 
modulus  function  was  determined  from  the  relaxation 
modulus  data  using  a  curve  fitting  technique.  Using  this 
relaxation  modulus  function,  the  effective  eleustic  modu¬ 
lus  wa^  calculated  by  carrying  out  the  linear  viscoelastic 
convolution  integral  for  a  strain  history  proportional  to 
the  motor  pressure  during  a  71  °C  ignition  pressurization 
with  a  36-ms  duration.  The  procedure  used  is  described 
in  detail  in  Reference  4.  The  resulting  modulus  values 
(accounting  for  the  error  estimate  produced  by  the  curve 
fit)  were 

+  3  Sigma  :  E  =  2.ZhMPa 
Mean  :  E  =  1.76MPa 

—  3  Sigma  :  E  =  1.32  MPa 

Equation  (1)  is  plotted  in  Figure  13  for  the  three  differ¬ 
ent  modulus  values.  As  shown  in  Figure  13,  the  analysis 
would  predict  unstable  deformations  for  the  mean  mod¬ 
ulus  value  if  burnback  were  ignored.  However,  for  the 
burnback  geometry,  the  analysis  predicts  only  a  small 
chance  of  instability  (i.e.,  the  modulus  must  be  at  or 
near  the  lower  3-sigma  value  to  cause  instability). 

CONCLUSIONS 

This  analysis  confirmed  the  potential  for  failure  of  the 


Sparrow  Mark  52  rocket  motor  through  large  bore  defor¬ 
mations  leading  to  choked  fiow  at  the  aft  end  of  the  pro¬ 
pellant  grain,  supporting  the  results  obtained  previously 
by  other  investigators.  The  following  specific  conclusions 
are  drawn  based  on  our  results: 

(1)  For  the  motor  operating  conditions  and  propellant 
properties  used  in  this  analysis,  the  structural  de¬ 
formation  is  marginally  stable.  Lower  propellant 
modulus  or  higher  burn  rate  (producing  a  higher 
head-end  pressure)  could  trigger  unstable  deforma¬ 
tion.  A  lower  liner  modulus  would  also  increase  the 
potential  for  instability. 

(2)  Because  of  the  significantly  higher  deformations,  the 
three-dimensional  structural  model  leads  to  a  sub¬ 
stantially  lower  margin  of  stability  than  would  be 
obtained  with  a  two-dimensional  structural  model. 

(3)  The  TEXGAP-2D/Approximate-3D  computer  code 
closely  approximated  the  radial  deformations  ob¬ 
tained  in  the  aft  end  of  the  motor  from  the  TEX- 
GAP-3D  code.  This  result  both  corroborates  the 
TEXGAP-3D  results  and  shows  the  usefulness  of 
the  Approxi-  mate-3D  feature  of  TEXGAP-2D  for 
analyses  involving  three-dimensional  geometry. 
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